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Using a relativistic mean-field single particle knock-out model for (e, e') reactions on nuclei, we 
investigate approximate treatments of Coulomb distortion effects and the extraction of longitudinal 
and transverse structure functions. We show that an effective momentum approximation (EMA) 
when coupled with a focusing factor provides a good description of the transverse contributions to 
the (e, e') cross sections for electron energies above 300 MeV on 208 Pb. This approximation is not 
as good for the longitudinal contributions even for incident electron eneriges above 1 GeV and if one 
requires very precise extraction of longitudinal and transverse structure functions in the quasielastic 
region it is necessary to utilize distortion factors based on a nuclear model and a more accurate 
inclusion of Coulomb distortion effects. 

PACS numbers: 25.30.Fj 



I. INTRODUCTION 



There continues to be considerable theoretical and experimental interest in extracting longitudinal and transverse 
structure functions as a function of energy loss for fixed three momentum transfer for a range of nuclei Q . For low 
to medium electron energies (200 MeV < E < 600 MeV) the (e,e') cross section is significantly affected by the static 
Coulomb field of the target nucleus for Z > 20 and it is necessary to use some method of removing the so-called 
"Coulomb distortion effects" in order to investigate the underlying physical processes in quasielastic scattering. For 
a precise extraction of the longitudinal response even at incident electron energies of 2 GeV, some correction for 
Coulomb distortion effects is needed. It would be very appealing if this removal of Coulomb effects could be effected 
by shifting the value of the experimental energy or scattering angle so as to compensate for the Coulomb distortion 
2] . While it is quite clear that this cannot be done exactly, the question is can one find some approximate treatment 
of Coulomb distortion effects that improves with increasing electron energy. 

In this paper, we investigate this question and report very good success in handling the Coulomb distortion effects 
in the transverse parts of the quasielastic cross section, but with less success in the longitudinal portion of the cross 
section unless we utilize a nuclear model of the transition current arising from the quasielastic knocking out of nucleons. 

In order to frame the discussion, it is useful to note that in the plane wave Born approximation (PWBA), the 
inclusive (e, e') cross section for electrons or positrons (assuming nuclear current conservation) is given by 

Mg,c.)j, (1) 

cos 2 — 

where q 2 = uj 2 — q 2 = — Q 2 is the four-momentum transfer, om is the Mott cross section given by <jm — (^fr) 2 4 e c , 

and Sl and St are the longitudinal and transverse structure functions which depend only on the momentum transfer 
q and the energy transfer lo. As is well known, by keeping the momentum and energy transfers fixed while varying the 
electron energy Ei and scattering angle 9 e , it is possible to extract the two structure functions with two measurements. 
The longitudinal and transverse structure functions in Eq. Q are squares of the Fourier transform of the components 
of the nuclear transition current density integrated over outgoing nucleon angles. Explicitly, the structure functions 
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for knocking out nucleons from a nucleus are given by 

S L (q,Lj) = Y, S ^PpJ \ N o\ 2d ^P (2) 

OtbS p 

S T {q,oj) = J2 Sa >Pp f (\N x f + \N y \ 2 )dn p (3) 

a b s p 

where the nucleon density of states p p = F^Xi > s p is the spin-projection of the continuum nucleon, the z-axis is taken 
to be along q, and S ab contains the spectroscopic and averaging factors for each bound orbital. The Fourier transform 
of the nuclear current J A '(r) is simply, 

= J J^(r)e iqr d 3 r. (4) 

and the continuity equation has been used to eliminate the z-component (N z ) via the equation N z — —^N . As 
noted above, the PWBA calculation breaks down for cases when the static Coulomb potential at the surface of the 
nucleus is not negligible when compared to the incident and outgoing electron energy. To set the scale, note that the 
Coulomb potential at the surface of 208 Pb is about 20 MeV. 

Based on our previous investigations of inclusive quasielastic scattering and a new appreciation of a focusing factor 
13- 0- 1 we investigate an improved effective momentum approximation with focusing (EMA-f) in this paper. 
Our goal is to seek a procedure for extracting the longitudinal and transverse structure functions from experimental 
data with minimal use of nuclear models. 



II. COULOMB DISTORTION EFFECTS IN INCLUSIVE QUASIELASTIC SCATTERING 

As discussed in a previous paper 0, we found a treatment of Coulomb distortion, labelled approximate DW, for 
(e, e') from medium and heavy nuclei that agrees to within a few percent to a full DWBA partial wave analysis 
which includes the static Coulomb field of the target exactly in the electron (or positron) wavefunctions. In this 
approximation we define an r-dependent momentum for the incoming and outgoing electron wavefunctions, *(r) = 

(pij — i J Q V c (r)dr)pij , where V c (r) is the potential energy of the lepton moving in the static Coulomb field of the 
target nucleus. This implies an r-dependent momentum transfer q'(r) = p^(r) — pj(r). The choice of an r-dependent 
momentum allows the distorted partial wave solutions of the Dirac equation with the static Coulomb potential present 
to be approximated by spherical Bessel functions with argument x = p'(r)r quite well from the origin out to more 
than 3 times the nuclear radius R. An additional requirement of approximating Coulomb distorted waves is to 
incorporate the Coulomb scattering phase shifts into the problem. As discussed in previous papers 0, we were 
able to achieve this by fitting the phase shifts as a function of the square of the Dirac quantum number k for the 
incoming and outgoing electron energies for the nucleus under investigation and then to replace k 2 by the classical 
angular momentum squared (r x p) 2 so that "plane- wave-like" lepton wavefunctions can be constructed that contains 
the effects of Coulomb distortion in the local momentum and the parametrized phase shifts which lead to what other 
authors refer to as focusing. Using this "plane-wave-like" wavefunction, we obtained an approximate "M0ller-like" 
potential given by, 

A a PP ro.DW (r) = j[8 i+ S,] fi i(A« - A,) ^ { r) ., - (5) 



where the phase shift parametrization is given by 
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where k 2 = (rxp) 2 and p is the electron momentum (for the initial and final kinematics) and we take the nuclear 
radius to be given by R = 1.12A 1 / 3 — 0.86A -1 / 3 . The two constants ao and 02 are fitted to two of the elastic scattering 
phase shifts (k = 1 and K = Int(pR) + 5) for the incident and final electron energy. The parameter A = a[p-f]L 2 
denotes a small higher order correction to the electron wave number which we have written in terms of the parameter 
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angular dependence of the vector r in the expression for the phase shifts. However, when combined with the nuclear 



a = —aZ( — p ' c ) 2 . The potential in Eq. JSJ) cannot be easily decomposed into a multipole expansion due to the 
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transition current density J M (r) as in Eq. Q, the modified Fourier transform can be obtained by three dimensional 
integration since the volume is limited to a sphere with radius of 3-4 times the nuclear radius R. We confirmed that 
use of the potential given in Eq. © reproduces the full DWBA results for the cross section very well for electron 
energies above 300 MeV and for momentum transfer greater than about 250 MeV. With this approximate £Wfour 
potential Au it is straightforward to calculate the exclusive (e,e'p) cross sections and modified structure functions. 
We showed 6] that using this new phase shift we can reproduce the full DWBA cross sections for (e, e'p) from medium 
and heavy nuclei very well. 

Again as noted in our previous paper [?J, this approximate DW potential which includes the local value of the 
potential (i.e, a function of r) and the focusing effect due to the phase shifts is very time consuming for the inclusive 
(e, e') reaction since we need to integrate over the solid angle of the outgoing nucleons. Thus, we proposed making 
further approximations for the inclusive process. In order to allow a straight-forward multipole expansion, we must 
remove any angular dependence from the phase shifts in Eq. JSJ. We investigated various methods of achieving this 
and found that neglecting the phase shifts entirely, but including a focusing factor was sufficient for the transverse 
contribution. However, this was not sufficient for the longitudinal term so we chose to include some of the effects of 
the phase shifts by averaging over the angles of the vector r in the phase shift parametrization. With these further 
approximations we were able to write the (e, e') cross sectionQ as 
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where the Fourier transforms of the transition current in S' L and S' T are replaced by 
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where < Sij > denotes an average of S(k 2 ) over the angles of the vector r. That is, the argument of the parametrization 
for the phase shifts in Eq. © is given by < n 2 f >=< (r x Pij) 



Jpi >/ )/4. In addition, we fix the 
direction of q' to be equal to the asymptotic momentum transfer, but we use q' 2 (r) = p' 2 (r) +p' 2 (r) — 2p\ (r)p'j: (r)cos9 e 
for the magnitude. We confirmed that for the kinematics under consideration this change is negligible. Using a toy 
model and the full three dimensional integration for the longitudinal terms and by comparing to the full DWBA for 
the cross section we found that this so-called ad-hoc model works very well for (e, e') on medium and heavy nuclei. 

However, use of an r-dependent momentum transfer requires a nuclear model to extract longitudinal and transverse 
structure functions. The challenge is to find some approximation with a constant shift of momentum transfer which 
approximates Coulomb distortion. One such approximation, referred to as the effective momentum approximation 
(EMA) replaces the r-dependent momentum for the incoming and outgoing lepton wavefunctions with a fixed value 
given by p' = p — V C (R C ) where R c is usually taken to be equal to or the nuclear radius R. Then one calculates the 
effective momentum transfer in terms of these momenta and the electron scattering angle. It is convenient to calculate 
V C (R C ) for a uniform charge distribution of radius R = 1.12A 2 / 3 — 0.86^4 -1 / 3 containing a charge of Ze. The result for 
R c < R is simply, V C (R C ) = — -f§r(3 — R 2 /R 2 ) for electrons. However, this approximation does not include the effect 
of the phase shifts which result in what other authors refer to as focusing. As discussed in some detail by the Basel 
group (3 , the effects of focusing can be included by multiplying the potential arising from the EMA wavefunction by 
the factor £^(0)/pi. We will refer to this approximation as EMA-f and we note that the (e, e') cross section for EMA-f 
is simply this factor squared times the EMA cross section. Thus, the EMA-f cross section is given by 
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The structure functions depend on the standard Fourier transforms of the nuclear transition current except that they 
are calculated as a function of q' rather than q. Note that Q' 2 = q' 2 — uj 2 . Since the sign of V c changes for positron 
induced reactions, clearly one can shift Ei for positron induced reactions as compared to electron induced reactions 
such that q' e _ — q' e+ and therefore the quantity ^u ) EMA ~ f /( a M(^p-) 2 ) is the same for electrons and positrons 
in this approximation. That is, the focusing factor must be removed before the overall structure functions are equal 
for electrons and positrons. 

We have investigated the validity of the EMA-f approximation for different choices of R Cl the argument of V c , by 
calculating the longitudinal and transverse contributions to the cross section using a relativistic nuclear model that 
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we have successfully used to describe a great deal of (e, e') data 0,0, El- In Fig. ^ we compare the cross section for 
the inclusive reaction (e, e') on 208 Pb for an incident electron energy of 310 MeV at a backward angle (6 — 143°). 
We show three curves in addition to diamonds representing the full DWBA calculation. As we have noted before, 
even at these low energies, our ad-hoc result agrees quite well with the full DWBA results. We show two EMA-f 
results. In one case we evaluated the Coulomb potential energy at the origin and in the other at two thirds of the 
nuclear radius which is in better agreement with the full DWBA results. Note that the cross section at such a large 
angle is dominated by the transverse contributions. In Fig. [21 we compare our ad-hoc results to EMA-f for the same 
two choices of the Coulomb potential for the longitudinal (upper panel) and transverse (lower panel) contributions 
to the cross section for 485 MeV electron on 208 Pb at a scattering angle of 60°. Note that the longitudinal and 
transverse contributions are of similar magnitude, but that the transverse cross section is much better described by 
the EMA-f than the longitudinal cross section. And furthermore, using the Coulomb potential at 2/3i? to calculate q' 
is a somewhat better approximation than using the origin value. In Fig. [3] we show the full cross section for this case 
as compared to the full DWBA calculation (where we cannot separate out the longitudinal and transverse terms) and 
we note that the ad-hoc result is very good and the EMA-f does not look so bad since over half the cross section comes 
from the transverse term which is reasonably well described by EMA-f. Of course if you are using EMA-f to make a 
Rosenbluth separation, the differential quality of the description of the longitudinal and transverse terms would lead 
to large errors. 

In Fig. 21 we examine the two separated contributions as in Fig. |3 except that we have increased the electron energy 
up to 800 MeV. Clearly the EMA-f approximation is much better at this higher energy although the discrepancy in 
the longitudinal case is still almost 10% at the quasielastic peak. In Fig. [SJ we increased the electron energy up to 2 
GeV while reducing the scattering angle to 20°. In all cases we use 208 Pb as the target. The EMA-f is almost exact 
for the transverse contribution to the cross section (lower panel) while the longitudinal contribution (upper panel) 
continues to have problems. 



III. PROPOSED SOLUTION AND CONCLUSIONS 



Our results show that even at rather high electron energies, the longitudinal contributions to the quasielastic cross 
section are not well described by the EMA approximation even with a focusing factor. Based on these results, we do 
not believe that a model-independent EMA-like approach can be used to extract the longitudinal structure function 
in quasielastic scattering. However, as the electron energy increases the EMA-f approach does get better and better 
so perhaps it can be used as the basic analysis tool, but with some model dependent corrections to the Coulomb 
distortion effects. In Fig. we show the ratio of the contributions to the quasielastic cross section for 485 MeV 
electrons on 208 Pb at a scattering angle of 60° calculated using our ad-hoc model over the EMA-f calculation using 
the Coulomb potential at |R for the longitudinal and for the transverse contributions. As was clear in Fig. [21 the 
distortion factor Dt for the transverse contributions differs from one by only a few percent over the quasielastic peak. 
However, the distortion factor Dl varies considerably from one at this energy. In Fig. [7|we repeat this calculation 
for incident electrons of 800 MeV. At this higher energy, Dt differs from one by less than 3% across the quasielastic 
peak while Dl deviates from one by up to about 7%. Based on these results and our results at higher energies, we 
propose that the quasielastic scattering cross section for (e, e') be written as, 
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where the distortion factors D^ and Dt are given by 
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Note that in the factor before the transverse structure function in Eq. IjlOjl . we are using the factor Q 2 /2q 2 rather 
than primed values so as to more easily define the distortion factor Dt- Clearly by using distortion factors Dl and Dt 
calculated from any reasonable nuclear model, Eq. (|10|) can be used in a Rosenbluth mode to extract the longitudinal 
and transverse functions as a function of q' and w from experimental data. 

In conclusion, we have shown that the effective momentum approximation (using the Coulomb potential at |R) 
with an overall focusing factor of (p^(0)/pi) 2 is a very good approximation of the Coulomb distortion effects for the 
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transverse contributions to the quasielastic cross section. However, for electron energies less that about 600 MeV it 
is not a good approximation of the longitudinal contributions. At higher electron energies this approximation does 
get better for the longitudinal contribution, but for accurate extraction of the longitudinal structure function it is 
necessary to use distortion factors calculated from a nuclear model. The procedure we have proposed minimizes 
the model dependence by only using the model to evaluate the Coulomb distortion effects not included in EMA-f. 
Based on the cases we have examined, we believe the errors in the extracted structure functions arising from Coulomb 
distortion effects using our proposed procedure should be less than 5% for incident electron energies above 600 MeV. 
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FIG. 1: Theoretical quasielastic scattering cross sections for 208 Pb with incident electron energy of Ei=310 MeV and electron 
scattering angle 6 — 143° as a function of energy transfer. The diamonds are the result of our full DWBA partial wave 
calculation. See text for details of the models and for details of the ad-hoc and EMA-f results. 
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FTO. 2: Theoretical onasielastic scatterine' nartial cross sections for 208 Pb with incident electron enerev of E,=485 MeV and 
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FIG. 3: Combined longitudinal and transverse contributions for the kinematics described in Fig. Q 



9 



4 



> 

0) 



in 



G 
3 



b 



3 - 



2 - 



> 

CD 



G 
3 



b 




200 

12 
10 

8 

6 

4 
2 



1 1 1 1 1 1 1 1 1 1 

E = 


i i i i 

800 MeV 


= 

< >-* . w 


60° 


— 

<•/ \ 

// \X 






7 \ 

y 

'/ 
// 


V. 

\C., 
V'- 


J* 11 IPs TUT" T~\ A 

/.' ad -hoc DWBA 


V* 
X* 

x* — 


EMA-f V(0) 




__. EMA-f V(2R/3) 

■ i i i 1 i i i i 1 





200 



300 



400 



500 



~~i 1 1 r~ 



—\ 1 1 1— 



| 1 1 1 r~ 

E=800 MeV 
= 60° 




ad-hoc DWBA 
EMA-f V(0) 
EMA-f V(2R/3) 



300 400 

cj (MeV) 



500 



FTO. 4: The miasielastie longitudinal and transverse contributions to the cross section for Pb with incident electron enerev 



50 
40 

30 

20 

10 

1 

90 
80 

70 

60 

50 

40 

30 

20 

10 
1 



E = 2.Q GeV 
= 20° 



ad-hoc DWB, 



EMA-f V(0) 
EMA-f V(2R/3) 



E = 2.Q GeV 
= 20° 



ad-hoc DWBA 



EMA-f V(0) 
EMA-f V(2R 



/3) 



200 



300 
cj (Me 



400 



Same as Fiff. Hlexcent that E,= 2.0 GeV and the electron scattering anrie 9 



11 




FTC1. 6: The ratio of the longitudinal and transverse (Dr and DtI eontribntions to the onasielastie eross section calculated 
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FTO. 7: The ratio of the longitudinal and transverse Dt contributions to the oua,sielastie cross section calculated usine the 



